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PREFACE 


This report was prepared in the course of a continuing study of 
Soviet research and development of high-current, high-energy charged 
particle beams and their scientific and technological applications. The 
report examines Soviet research on the propagation of intense relativis¬ 
tic electron beams (IREBs) through high-pressure air and gases 
(p > 10“^ Torr), as reported in Soviet open-source technical publica¬ 
tions. A companion RAND report, Soviet Research on the Transport 
of Intense Relativistic Electron Beams Through Low-Pressure Air, 
R-3309-ARPA, August 1986, by the same author considers Soviet 
research on the propagation of IREBs through low-pressure air and 
gases (10~® Torr < p < 10~® Torr). 

The report should be of interest to specialists in pulsed-power and 
high-current charged-particle-beam research. 

This study was sponsored by the Defense Advanced Research Pro¬ 
jects Agency, Office of the Secretary of Defense (OSD), and was carried 
out in the Applied Science and Technology Program of the National 
Defense Research Institute, RAND’s OSD-supported Federally Funded 
Research and Development Center. 
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SUMMARY 


The Soviets have actively pursued research into the propagation of 
intense relativistic electron beams (IREB) (beam energies above 450 
keV and beam currents in the tens of kA range) through air and other 
gases since the early 19708. The present report analyzes the Soviet 
developments of IREB propagation through background air at pres¬ 
sures from 10'^ Torr to atmospheric, as reflected by Soviet open source 
technical literature published over the past 15 years. A previous report 
by the same author considered Soviet developments on IREB propaga¬ 
tion through air at low pressures (10^® to 10“^^ Torr) [1]. 

Soviet research on IREB propagation in air has greatly increased 
recently, as evidenced by: 

• A large increase in the number of researchers and published 
papers on IREB propagation from 1980 to 1984. 

• A coincident disappearance from publication in any area of sci¬ 
ence of five important experimental research teams, suggesting 
censorship action imposed on these teams. 

The most important findings of Soviet research seem to match U.S. 
observations regarding propagation, as follows: 

1. The formation of a plasma channel created by an IREB propa¬ 
gating through background air and the effect of beam parame¬ 
ters upon the plasma channel parameters (and vice versa). 

2. Determination of the background air pressure for the optimum 
transport of IREB. Two optimum pressure ranges were deter¬ 
mined: 

(a) The “ion focused regime”—0.06 to 0.09 Torr 

(b) The “low pressure window”—1 Torr 

3. Observation of current enhancement, whereby the IREB-induced 
current in plasma is higher than the initial beam current. 

4. The effect of resistive hose instability upon IREB propagation. 

Both the United States and the USSR are concentrating on the 
same effects inherent to IREB propagation through air. However, an 
important difference between the Soviet research, as presented in the 
open literature, and the U.S. research is that the Soviet high current 
experiments are carried out with pulsed diode machines in the 1 to 1.5 
MeV range and with a poorer quality beam than the U.S. experiments 
such as on the ATA 50 MeV accelerator and other beams. There is a 
striking absence of any high energy, high current research in the 
USSR. 
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Most Soviet research has been carried out at nine institutes, five of 
which have made extensive experiments using IREBs with energies of 
0.5 to 1.5 MeV, beam currents of 15 to 50 kA, and pulse lengths of 20 
to 60 ns, with beam propagation lengths from fractions of a meter to 3 
meters. In one case, experiments using beam energies of 2 to 20 MeV 
were made at low beam currents. 

Soviet theoretical investigations on the formation of a plasma chan¬ 
nel by the passage of an IREB through air have been carried out by 
computer simulation made on the BESM-6 computer for cases of a 1 
MeV, 30 kA IREB. A new model used was more adapted to higher 
pressures and longer distances. This model took into account changes 
in plasma channel parameters due to IREB and changes in beam 
charge and current density profiles. Pressures considered in this model 
were from 1 to 100 Torr and the electron energy inside the plasma 
channel was tens of eV. 

According to Soviet researchers, for optimum beam transport 
through air it is imperative to ensure conditions that allow full neutral¬ 
ization of the IREB^s self-fields along the entire path of the beam’s 
transport. The propagation is determined by the plasma channel con¬ 
figuration while the trajectory is determined by movement of electrons 
that have lower energy than the main component of the electron beam. 

Soviet emphasis on IREB propagation through fairly high pressure 
may be important in light of the recent increased effort. The high 
pressure range considered in this r^ort is basically applicable to iner¬ 
tial confinement fusion, where the beam must be transported with 
minimum losses from the accelerator to the target. However, in the 
case of IGF, the gas through which the beam propagates usually is not 
air, because propagation throxigh air rather than through other gases 
(especially He) was found to be much less efficient. This pressure 
range also applies to the passage of IREB through the atmosphere and 
is pertinent to endoatmospheric charged-particle beam weapons and 
other military applications. 
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GLOSSARY 


The IREB and plasma channel parameters being considered by the 
Soviet researchers are as follows: 


IREB Parameters 

Symbol 

Units 

Beam energy 

Eb 

MeV, keV 

Beam current 

h 

kA 

Beam current density 

Jb 

kA/cm^ 

Beam radius 

n 

cm 

Beam radius (at input to drift tube) 

rbo 

cm 

Azimuthal magnetic field 

Be 

G, kG 

Time after pulse 

t 

ns, iis 

Plasma Channel Parameters 

Background gas pressure 

P 

Torr 

Optimum background gas pressure 

Popt 

Torr 

(for best IREB propagation) 

Propagation length 

Z 

m 

Total length of drift tube 

L 

m 

Gas or channel plasma electron density 

rie 

cm~^ 

Gas or channel plasma ion density 

Hi 

cm~^ 

Characteristic plasma electron energy 


eV 

Conductivity of plasma channel 

a 

mho/m 

Plasma current 

Ip 

kA 

Plasma current density 

Jp 

kA/cm^ 

Reverse plasma current 

Ir 

kA 

Net total current 

It 

kA 

Drift tube inside radius 

Rc 

cm 

Plasma channel radius 

^P 

cm 

Plasma magnetic field decay time 


ns 


xvii 



I. INTRODUCTION 


Soviet research on the propagation of intense relativistic electron 
beams (IREB) through fairly high-pressure air (pressure range 0.1 to 
760 Torr) since the early 1970s included the study of the plasma chan¬ 
nel created by the passage of the electron beam through air, the resis¬ 
tive hose instability and its effect on beam propagation, the effect of 
self-fields, current enhancement, gas expansion, return currents, 
inherent beam energy spread, and other factors. 

This report covers Soviet developments in IREB propagation 
through air where the beam is not focused by external magnetic fields. 
The information was obtained from Soviet open publications with 
emphasis given to the last ten years of beam propagation in the USSR. 
The background pressure of air considered here varies from 0.1 Torr to 
atmospheric pressure. Soviet work on IREB transport through air at 
low pressures (from 10“® to 10“^ Torr) was considered in a previous 
report [1]. Earlier work on Soviet developments of IREB propagation 
through air can be found in an earlier RAND publication by Kassel 
and Hendricks [2]. 

Most Soviet beam propagation research has been carried out at nine 
institutes as shown in Table 1. Researchers and their leaders are listed 
under their institutes in the appendix. The chronology of the key 
Soviet research on IREB propagation through air is presented in Fig, 1 
for the nine institutes and their leading researchers. 

Figure 1 also indicates the extent, initiation, and termination of the 
specific research within the individual institutes. Note the short time 
span of the experimental research carried out by the Yu. P. Usov team 
at lYaF-TPI where many researchers apparently made a concentrated 
effort, and the similar short and limited effort made at NIIEFA and 
NIITP but with no recent publications. Also of interest is the large 
effort carried out by two teams at KhFTI, with a complete cutoff in 
publications in 1981 and 1982, respectively. With the termination of 
publication on these five important experimental investigations comes 
the onset of five research efforts that are currently represented by 
recent publications in lYaF-TPI (Team 1: V. P. Grigor’yev), lAE (L. I. 
Rudakov), FIAN-IOF (A. A. Rukhadze), VEI (A. V. Zharinov, M. A. 
Vlasov), and MFTI (K. V. Khodatayev). 

Seven of the abovementioned institutes carried out experimental 
investigations on IREB propagation through air. Table 2 shows the 
accelerators used in this research, together with the beam and propaga¬ 
tion parameters. 
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The experimental research carried out in the first five institutes 
shown in Table 2 utilized IREBs with energies of 0.5 to 1.5 MeV, beam 
currents of 15 to 50 kA, and pulse lengths of 20 to 60 ns, with beam 
propagation lengths from fractions of a meter to three meters. In the 
case of the KhFTI research, two experiments are included where higher 
beam energies were used but at low beam currents. 

Theoretical work of any consequence was carried out only at lYaF- 
TPI, MFTI, and in a very limited way at VEL Out of the seven insti¬ 
tutes that have performed experiments, only lYaF-TPI has carried out 
computer simulation work of any magnitude. The theoretical work at 
MFTI is somewhat of a mystery, the papers not being clear on applica¬ 
tions and analysis. However, the researchers at MFTI under 
Khodatayev’s leadership are of interest, and their citations refer to 
many U.S. researchers working on IREB propagation. 

Figure 2 (a and b) plots the total number of Soviet open source 
papers published specifically on research of IREB propagation through 
and injection into air for the past 15 years, and the number of authors 
of these papers. The research papers on the use of IREBs for optical 
pumping, beam welding, and other electron beam applications that are 
not obvious for beam transport are not included. The number of 
papers as a function of years gradually increases with distinct peaks 
during 1974, 1978, 1981, and 1984. The great jump of activity in 1984 
demonstrates an increase of more than a factor of three of papers pub¬ 
lished in 1984 over the number published in 1983. 

The plot of number of authors publishing the research papers is a 
demonstration of relative effort. Each author was considered to be 
active (and was therefore counted) a year before the date of his first 
published paper until one year after his last published paper. Figure 2b 
demonstrates a sharp increase in the number of authors involved in the 
IREB transport research from 1980 to 1984. 

This abrupt jump in the number of publications and authors partici¬ 
pating in the research in 1984 may signify an increase in the recent 
overall Soviet effort on IREB transport through air, but it may also be 
the result of a Soviet declassification of earlier research, which is only 
now able to be published. The peak may therefore represent the 
“dumping” of papers. 

The existence of an active Soviet CPBW program cannot be directly 
inferred from the available data as presented in the Soviet open source 
literature because of the low energies (of 1.5 MeV and below) and beam 
propagation distances (of 3 meters or less) used in the experiments. 
The Soviet work, however, does demonstrate a broad research effort, 
involving some key elements of IREB propagation through the air. 
The sharp increase in published Soviet papers in recent years implies a 
significant increase in the research effort in this area. 
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Years 

Fig, 2—Number of Soviet open source papers (lower curve) published 
on IREB propagation through air, and the number of authors of 
these papers (upper curve) as a function of years 


Number of papers 


II, lYAF-TPI RESEARCH 


THEORY AND COMPUTER SIMULATION 

Propagation research of the Nuclear Physics Institute attached to 
the Tomsk Polytechnical Institute (lYaT-TPI) has been performed by 
two separate research teams. The first, headed by V. P. Grigor’yev, 
carried out theoretical calculations and computer simulation studies of 
IREB propagation from 1977 to the present. The second team, headed 
by Yu. P. Usov, performed several experiments with IREB propagation, 
described in the next subsection. 

In 1979“1980 the theoretical team carried out computer simulation 
studies using the BESM-6 computer to investigate plasma channels 
formed by the injection of IREBs into air and neutral gas at pressures 
in the range of 1 to 100 Torr. The team considered plasma channel 
inhomogeneities along the beam axis, the effect of the Bz field gradient 
along the axis and its effect upon beam radius, and beam current neu¬ 
tralization in the plasma channel for beam pulse lengths of 10 to 80 
nsec. For these calculations the researchers took the values of the 
beam parameters from early transport experiments using the “Tonus” 
accelerator with a beam energy of 1 MeV and beam currents of 20 to 
30 kA [3]. 

The computer simulation studies as carried out in 1979-1980 and 
described in Refs. 4 and 5 did not take into account the changes in 
charges and electron beam current distribution within the beam during 
IREB transport. The data and its discussion as presented in Figs. 3 
through 12 are based on the simpler models, which are applicable only 
at low background air pressures and short beam transport distances 
and have been investigated in the pressure range from 0.2 to 40 Torr. 
Computer simulation investigations carried out in 1983 as described in 
Ref. 6 presented a more sophisticated model taking into account the 
beam’s charge and current distribution. These more recent data are 
presented in Figs. 13 through 20 and were determined in the back¬ 
ground air pressure range from 1 to 100 Torr, 

The parameters involved in the above computer simulations were: 
the beam current, 7^; total system current /y [7r = 7^ -f- 7p, where Ip 
is the plasma current]; beam current density, the azimuthal mag¬ 
netic field Be, beam radius, electron plasma density, characteris¬ 
tic energy of plasma electrons and the conductivity of the plasma 
channel, a, as functions of background gas pressure, p; propagation 
length, z] time after the pulse, t; beam radius, r^; and drift tube inside 
radius, Rc. (See the Glossary for symbol definitions.) 


7 



8 


The variation of Be and It along the beam axis for background air 
pressures of 1 and 10 Torr are shown in Fig. 3, and the axial variation 
of as a function of p is shown in Fig. 4. The radial Be field distribu-- 
tion for a background pressure of 1 Torr is shown in Fig. 5 [4], 



Fig, 3—Axial profile of the azimuthal magnetic field B^(solid lines) and 
the total current Irji (dashed lines) for p = 1 and 10 Torr for times t = 
50 ns - curve 1, 30 ns - curve 2, and 10 ns - curve 3 [4] 



Fig, 4—Axial profile of the azimuthal magnetic field in a 1 meter 
long drift tube, t = 30 ns [4] 


, kA 
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1 - .08 m 

2- 1 m 

3- 1.92 m 


Fig. 5—Radial field distribution for p = 1 Torr, t = 20 ns as a 

function of various axial distances [4] 


Soviet theory points out that the decrease in indicates a notice¬ 
ably defocusing condition of the beam at t > 30 ns; but at short beam 
transport distances, the basic structure of the electron beam does not 
change even at the high air pressure of 40 Torr. 

As seen in Fig. 5, the magnetic field intensity is located in a narrow 
band at the edge of the beam with the area in the middle of the beam 
practicaUy magnetically neutralized. This Be field radial distribution is 
the same for other background pressures, except that the degree of 
magnetic neutralization in the center of the beam decreases with an 
increase in pressure (as seen in Fig. 4). 

The longitudinal profiles of the plasma density in the plasma chan¬ 
nel at different pressures are shown in Fig. 6. The azimuthal magnetic 



Fig. 6—Plasma density longitudinal profile at p = 1 Torr (dashed 
line) and 10 Torr (solid line) at times of: t = 50 (1), 30 (2), 
and 10 ns (3) [4] 
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field Be decreases along the beam path for IREB tran^ort over dis¬ 
tances greater than 1 meter in the pressure range > 10 Torr. Such a 
Be gradient demonstrates a significant erosion of the beam front. This 
effect, however, can find an application in shortening the beam pulse. 
The simulation also showed that the radius of the plasma channel can 
be considerably larger than the beam radius at pressures below 1 Torr. 
At higher pressures the radii of the beam and plasma channel are prac¬ 
tically the same. The beam current density radial distribution as a 
function of pressure is shown in Fig. 7 [5]. 



1 — p« 1.10 Torr 

2 — p = 0.2, 0.5 Torr 

Fig. 7—Radial distribution of the beam current density for several 
background air pressures [5] 

The measured and calculated values of the total beam current pulses 
are shown in Fig. 8 for a z = 30cm [5], Experimental values of the 
total current are shown by dashed lines demonstrating good agreement 
between experiment and theory for these cases. 

The time dependence of the plasma density created by the transport 
of the electron beam through a 12 cm radius drift tube for different 
pressures is shown in Fig. 9 [6]. The plasma density produced from 
the collision and cumulative ionization in the area of the beam (r = 0) 
increases rapidly and reaches a maximum value of about 3 x 10^^ cm"’^ 
at a pressure of 10 Torr and 2 x 10^^ cm~^ at a pressure of 0.2 Torr. 

Outside the electron beam the plasma density increases only by 
cumulative ionization processes and is strongly dependent upon the 
background gas pressure increasing by more than two tens for a 
decrease in pressure from 10 to 0.2 Torr. The plasma density within 
the beam in the pressure range of from 0.2 to 10 Torr was weakly 
dependent upon pressure and mainly dependent upon beam electron 
density. 
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0 AO 80 
t, ns 


p = ,2 Torr 

1 - Rg = 4.6 cm 

2 - 12 cm 


p - 1 Torr 

1 - R^ = 4.6 cm 

2 - R^ - 6 cm 

3 - - 12 cm 


p - 10 Torr 

1 - R^j = 4.6 cm 

2 - R^ = 6 cm 

3 - R^ - 12 cm 

-* beam current 
dashed line - experiments 
with R^j - 4.6 cm 


Fig^ g—Total beam current pulse as a function of background pressure 
and drift tube radius [5] 



Fig. 9—Time dependence of the plasma density (n^) as a function of 
background pressure for r - 0 and r = 5 cm [5] 


12 


The radial distribution of the plasma density as a function of several 
pressures and several beam drift tube radii is shown in Fig. 10. In the 
0.2 to 10 Torr pressure range, the plasma density follows the radial 



= 4.6 cm, 1 - p = 0.2 Torr, 2 - 0,5 Torr, 3 - 1 Torr 
R^, = 6 cm, 4 - p = 0.2 Torr, 5 - 0.5 Torr, 6-1 Ton 
Rq = 12 cm, 7 - p = 0.2 Torr. 8 - 0.5 Ton. 9-1 Torr 

Fig. 10—Radial plasma density distribution (n^) as a function of 
background pressure and drift tube size for a time of 37 ns 
after initiation of the beam pulse [5] 


profile of the beam density, depends weakly upon pressure, and falls off 
rapidly beyond the beam with increased background pressure. The 
radial plasma density distribution was flat at the beam axis and steep 
at the beam edge. The dependence of the plasma channel radius upon 
pressure for several drift tube r adii i s shown in Fig. 11. 

The characteristic plasma electron energy (c^) averaged over the 
beam cross section, as a function of time for various background pres¬ 
sures, is represented in Fig. 12a. The radial distribution of the plasma 
electron energy is shown in Fig. 12b for the various background pres¬ 
sures. Figure 12a shows that c/g has a maximum at the start of the 
beam current pulse when the fields produced at the beam front are 
large and then it decreases with time. In the radial distribution of 
(Fig. 12b) the value of in the area of the electron beam remains 
fairly constant and falls off slowly with increasing radial distance. The 
value of decreases with increasing background pressure, which is due 
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A Ton 

Fig. 11—Plasma channel radius r as a function of background pressure 
for several drift tube radii [5] 



Fig. 12—Characteristic average plasma electron energy as a 
function of time (a) and its radial distribution (b) at time 
t = 37 ns and drift tube radius of 12 cm for various 
background pressures [5] 
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to the increase of the collision frequency i/gm the plasma electrons 
with the gas molecules as the pressure increases. The plasma electrons 
lose their energy with each collision in time 1/vem and do not have 
enough time to obtain a high energy. 

The beam was found to induce a current flowing opposite to the 
direction of the electron beam in the plasma. Most of this current 
flows within the electron beam; and its current density, Jp, was found 
to decrease rapidly with increasing radius. 

The lYaF-TPI researchers state that they are quite satisfied with 
the mathematical model as presented in Refs. 5 and 7 that demon¬ 
strated good agreement between experiment and theory in the pressure 
range of 0.2 to 10 Torr. The investigators determined an optimum 
beam neutralization at the pressure of about 1 Torr and found that it 
is very dependent upon the inner radius of the metal tube. They found 
that this optimum beam neutralization increases with an increase in 
drift tube radius. 

The theoretical models of Refs. 4 and 5 allowed the lYaF-TPI 
researchers to narrow the gap between theory and experiment. How¬ 
ever, in these models the beam characteristics were predetermined and 
the change in charge and current distributions within the beam were 
not taken into account in the process of beam transport. These 
assumptions, as was shown in the experiments of Refs. 3 and 7, are 
correct only in the case of low pressures and short distances, because at 
higher background gas pressures and larger beam transport distances, 
the beam current and charge density profiles change, which leads, in 
turn, to a change in the plasma channel parameters. Another model 
for use at higher pressure and longer transport distances was proposed 
in 1983 in Ref. 6, taking into account the change in the parameters of 
the plasma channel formed by the IREB and the changes in beam 
charge and current density profiles. 

A numerical simulation, made in Ref. 6, assumed an IREB through 
air at pressures from 1 to 100 Torr, with the electron energy in tens of 
eV inside the plasma channel, and approximating the recombination 
coefficient by the Mehr and Biondi formula [8]. The plasma density 
Tie as a function of transport distance and time is shown in Figs. 13 
and 14, respectively [6]. 

The plasma density was seen to reach 4.5 x 10^® cm~^ in 30 ns at 
pressures 50 to 70 Torr and about 2.5 x 10^^ at pressures from 1 
to 20 Torr. Beyond 30 ns the plasma density decreased because the 
recombination rate (mostly dissociative) was greater than the ioniza¬ 
tion rate. The plasma channel radius at pressures greater than 1 Torr 
coincided with the beam radius, and the radial profile n«(r, 2 ,t) was 
very close to The plasma channel was inhomogeneous along 
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ntlr‘0)^10'^\ CM* 



Dashed lines represent a pressure of 1 Torr 

1 - t = 10 ns 2 - t = 40 ns 

Solid lines represent a pressure of 20 Torr 

1 - t = 6 ns 2 - t = 10 ns 

3 - t = 17 ns 4 - t = 40 ns 

where r^^^ is the beam radius at the entrance to the drift 

tube and is the drift tube radius 


Fig. 13—The plasma electron density (n^) distribution along the 
drift tube with rj^^^ = 1 cm, = 2 cm [6] 

the length of the drift tube, especially at the higher pressures and for 
large initial beam densities. 

The dependence of the density of the plasma in the channel upon 
gas pressure and initial beam current density was found to be complex, 
because of all these factors affecting beam density and ionization rate. 
At short distances from the point of beam injection the plasma electron 
density /ig was greater at higher gas pressures. The dependence was 
different at the output of the drift tube, however. For t < 10 to 12 ns, 
Mg was lower at p == 1 Torr than at p = 10 to 100 Torr. However, for 
t > 12 ns, /ig increased more rapidly at p = 1 Torr. 

The characteristic energy of the plasma electrons (e^) was found to 
depend only weakly upon the gas pressure. The value of was found 
to vary from 1 to 14 eV over this pressure range and along the whole 
length of the channel. The higher values of 6^ (5 to 14 eV) were found 
for p = 1 Torr and the lower values (0.5 to 2 eV) for p = 100 Torr. 
The value of c* was found to be very weakly dependent upon the 
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Solid lines represent Z = 8 cm, with pressures of; 

1 - 1 Ton* 2-10 Ton* 3-20Torr 
4-50 Torr 5-70 Torr 

with r^j^ = 1 cm, = 6 cm 
6-20 Torr, rj^ = 3 cm, = 6 cm 
Dashed lines represent Z = 92 cm, with pressures of; 
1 - 1 Torr 2-10 Torr 3 - 20 to 70 Torr 
with r^j^ = 1 cm, R^ = 2 cm 
4-20 Torr, = 3 cm, R^ = 6 cm 


Fig. 14—Dependence of n^ upon time [6] 


radius. This can also be seen in Fig. 12 with earlier calculations of ek 
as a function of r and t [5]. 

The plasma channel conductivity (a) was found to have a strong 
dependence upon gas pressure and at pressures greater than 1 Torr was 
determined totally by plasma electron collisions with neutrals. This <r 
depends not only upon ek but also upon the ratio of the densities of the 
plasma and neutral molecules. The conductivity was thus found to 
decrease rapidly with an increase in gas pressure, and its radial distri¬ 
bution is similar to that of the rie distribution. The dependence of ek 
and <T upon gas pressure is demonstrated in Fig. 15. 

The simulation of the total current It ^ h Ip and the azimuthal 
magnetic field Be demonstrates that the basic plasma current flows in 
the area of the axis of the beam and decreases as r approaches 
With increase in gas pressure and because of the decrease in the 
plasma channel conductivity, the current neutralization weakens and 
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cr, mho/m £|^, eV 



Solid curves represent a and dashed curves represent 

1 - t = 10 ns, z = 8 cm 

2 -1 » 33 ns, z = 8 cm 

3 - t = 10 ns, z = 92 cm 

4 - t » 33 ns, z = 92 cm 


Fig. 15—Dependence of the characteristic energy of the plasma 
electrons (6|^) and the conductivity of the plasma channel (a) 
\: 5 )on gas pressure (p) for r = 0, Tbo = 1 cm, = 2 cm [6] 


Be grows near the point of beam injection. The radial distribution of 
current neutralization becomes more homogeneous. The Be radial dis¬ 
tribution found in Ref. 6 is similar to the earlier findings made in Refs. 
4 and 5. According to Ref. % Be increases with time, and the sharp 
dependence of as a function of z is seen only at the beginning of the 
drift tube. At the end of the tube the dependence on t becomes insig¬ 
nificant, as demonstrated in Fig. 16. 

In the light of the above theoretical results, the lYaF-TPI research¬ 
ers explain electron beam focusing as follows: A drift of plasma elec¬ 
trons is observed along the beam radius as a result of the movement of 
the electron beam in the Be, which increases in time and falls off with 
distance. Because Be increases with time, the electrons obtain a veloc¬ 
ity directed toward the beam axis, and because of the Be gradient along 
2 , the velocity is directed along the axis. Beam electrons follow trajec¬ 
tories that are oscillations with either decreasing or increasing radial 
amplitudes, depending upon which forces of the magnetic field predom¬ 
inate. Thus the beam either focuses or defocuses. 
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Dashed curves: p = 1 Torr 
1 - t = 7 ns 
2-t* 17 ns 
3 - t = 33 ns 
Solid curves: p = 50 Torr 

1 - t = 6 ns 

2 - t = 10 ns 

3 - t = 17 ns 

4 - t = 33 ns 


Fig. 16—Axial profile of the azimuthal magnetic field for r = = 

1 cm, Rj. = 2 cm [6] 


The beam dynamics calculations are shown in Figs. 17-19 [6]. In 
Fig. 17, curve 1 represents the beam pulse without any current loss at 
the output with the following two conditions: 

p = 1 Torr, r^o = 1 cm, Rc 2 cm 

p = 20 Torr, r^o = 3 cm, Re “ 6 cm 

Figure 17 indicates that a beam having a high current density of about 
7.6 kA/cm^ at the input to the drift tube can be transported without 
losses only for pressures less than 5 Torr. At higher pressures, only 
the leading edge of the beam pulse can pass, while the trailing edge is 
cut off. This shortens the pulse and decreases its amplitude at the 
drift tube output. The cutoff effect is more evident at higher gas pres- I 
sures and beam current densities, in longer drift tubes, and with * 
smaller tube radii. As an example, beam transport was computed 
for beam pulse increasing from = 0) = 23.7 kA, = 

7.6 kA/cm^ to Ii, max (2=0)- 50 kA, J^o = 16 kA/cm^. Pulse length 
was found to decrease from 33 to 26 ns for L = 1 m, J?c ^ 2 cm. The 
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Ib/Ib, 



1 - at the Input to the drift tube 

2 - 7 at the output of the drift tube 

2 - p = 5 Torr 

3 - p = 10 Torr 

4 - p = 20 Torr 

5 - p = 30 Torr 

6 - p = 50 Torr 

7 - p = 70 to 100 Torr 

with lb max = 23.7 kA. r^o = 1 cm. = 2 cm 

Fig. 17—Electron beam current pulse shape [6] 


ratios of the beam current values at the input and output of the drift 
tube were 0.56 and 0.38, respectively. 

The above theoretical results were in good agreement with experi¬ 
ments performed in 1967 [9] when it was determined that current was 
cut off at the trailing edge of the pulse in drift tubes 1 meter long at 
pressures above 6 Torr. 

Figure 17 indicates that the current cutoff stabilized at pressures 
from 50 to 70 Torr, and there was no further pulse shortening beyond 
70 Torr. There also was no current cutoff at pressures below 5 Torr. 
Both effects were attributed to current neutralization. At low pres¬ 
sures, the gradient of the axial magnetic field is flattened out (see Fig. 
16); and therefore there is no pulse current cutoff present. The con¬ 
ductivity of the plasma channel is then high enough (Fig. 15) and the 
level of current neutralization is also sufficiently high during the entire 
pulse. As pressure increased beyond 5 Torr, plasma channel conduc- 
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(a) p = 1 Torn 1 -1 ■ 10 ns, 2 -1 = 33 ns, 

3 - t = 65 ns 

(b) solid curves - p = 50 Torr: 1 - t = 10 ns, 

2 - t = 17 ns, 3 - t * 33 ns 
dashed curves - p = 70 to 100 Torr: 

1 - t = 10 ns, 2 - t = 33 ns 

(c) p = 20 Torr: 1 - t = 10 ns, 2 -1 = 17 ns, 

3 - t = 33 ns, 4 - t = 65 ns 


Fig. 18—Dependence of the beam radius upon the axial length z [6] 


tivity decreased (Fig. 15); and the level of current neutralization was no 
longer sufficient to flatten out the Bq field gradient and to ensure the 
beam transport without loss. Above 50 Torr, the conductivity of the 
channel decreased to such an extent that plasma current had practi¬ 
cally no effect upon the axial magnetic field distribution. Bq was then 
determined only by the inhomogeneity of the injected beam at the 
beam front and by its current density change in the process of beam 
transport. 

The change in beam parameters along the drift tube is shown in 
Figs. 18 and 19. At short distances from the beam injection plane 
(2 ~ 8 to 16 cm), the beam was focused by a sharp increase of Bq in 
time. As the beam traveled further, it started spreading; and the basic 
electron losses occurred at 2 ^ L/2. At the end of the drift tube, at 
pressures greater than 10 Torr and T > 18 ns, the electron beam den- 
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Rjj = 2 cm, p = 50 Toit_ 

1 - t = 10 ns, 2 - t = 14 ns, 3 - t = 17 ns. 
4 - t = 20 ns, 5 - t * 26 ns, 6 - t = 33 ns 


Fig. 19—Relative beam current losses for = 1 cm [6] 


sity on axis sharply dropped, and the radial beam profile rttir) was 
homogeneous along the entire cross section of the drift tube. Beam 
spreading decreased with increasing drift tube radius Rc (Fig. 20). This 
can be attributed to the increase of the inductance of the system with 
increasing Rc and the increase in current neutralization. 





1 - Rg - 2 cm, 2 - * 4 cm, 3 - = 10 cm 

, dashed line - R^, = 20 cm 

Fig, 20—Beam radius r^ as a function of distance and drift 
tube radius [6] 
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Soviet researchers stress that even dense electron beams have little 
beam loss at short distances from the beam injection plane up to pres¬ 
sures above 50 Torr (Fig. 19). A detailed analysis of the propagation of 
an IREB along a plasma channel with a return current (current flow¬ 
ing opposite to the IREB in the channel) was made in Ref. 10. Here 
the radial distribution of the charge density and currents of the elec¬ 
tron beam, plasma electrons, and ions have been determined for the 
beam radius smaller than, equal to, and larger than the radius of the 
plasma channel. 

The beam is assumed to be in equilibrium, with its space charge par¬ 
tially neutralized by the background gas ions. This steady state condi¬ 
tion is possible when an IREB in injected into low-density plasma 
(where nt is the density of the electron beam and n/ is the density of 
ions of the background gas). However, when the electron beam is 
injected into a dense plasma, space charge neutralization is also accom¬ 
panied by a partial beam current neutralization, characterized by 
plasma electrons flowing opposite to the direction of the electron beam. 

The duration of the return plasma current flow was found to be 0.1 
microseconds for 1cm, rig -- 10^"^ to 10^^, and electron plasma tem¬ 
perature of ~ 1 eV and is longer than the usual beam current pulse 
length [10], As a result, when Soviet researchers considered the equi¬ 
librium condition of electron beam transport, they found it necessary 
to take into account this plasma electron current as well as the effect 
of the background gas ions. 

The radial potential distribution as a function of the system configu¬ 
ration, the electron beam radius, the plasma radius, and the total beam 
current have been investigated using the equations developed in Ref. 
10. A specific calculation was made for the case when the electron and 
ion plasma radii were the same (n = = Tp) and when the IREB 

radius was equal to the plasma radius (rt = Tp). 

With a net negative charge excess along the beam axis, the beam 
and plasma electron density increased to the outer beam periphery 
while the ion density decreased [10]. However,. with a net positive 
charge excess along the axis, the beam and plasma electron density 
decreased to the beam periphery, while ion density increased. When 
the electron beam current density on the axis was greater than the 
plasma electron current density, the electron beam current density 
increased with an increase in beam radius while the plasma current 
density decreased. Thus, the radial distribution of the total current 
density of the system represented a hollow beam (Fig. 21a), which was 
in agreement with the 1971 experiments with IREB injection into 
dense plasma [11]. The current transported by the ions could be 
neglected in that case. 
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(a) (b) (c) 



a - = Tp, b - < Tp, c - Tj, > 




Fig. 21 “Radial distribution of the current density of the beam, plasma, 
and total current flow [10] 


When the IREB radius was smaller than the plasma radius 
(fp > Tf,), the current density of plasma electrons dropped to zero; with 
increasing radius and changing direction, it increased monotonically as 
it approached the plasma channel boundary (Fig. 21b). The radial pro¬ 
file of the total current in this configuration formed a double wall tube. 

A beam radius can be larger than a plasma radius (r^, > Tp) when the 
beam is transported along thin metal wires or along a thin laser beam. 
The electron beam current density was found in this case to increase 
slightly to the plasma boundary; beyond it increased more rapidly to 
the beam boundary as shown in Fig. 21c. 

The Soviet researchers state that IREB propagation in plasma with 
return current permits effective beam transport with beam currents 
that are considerably higher than those in static cases of neutralized 
electron beams. These beams have an application where it is necessary 
to create IREB with very high currents [10]. 

lYAF-TPI RESEARCH—EXPERIMENTS ON CHANNEL 
FORMATION AND HOSE INSTABILITIES 

The lYaF-TPI experimental research team headed by Yu. P. Usov 
performed several experiments in 1975-1977 involving plasma channels 
induced by IREB, hose instability, the optimal background pressure for 
IREB propagation, and the energy spread of beam electrons. These 
experiments were carried out using the Tonus accelerator operating at 
1 MeV with beam currents from 35 to 40 kA, and a pulse length of 60 
ns. No publications by this team were found after 1977. 

The formation of the plasma channel was first observed by Soviet 
researchers in 1972 when an IREB was injected into air and gases at 
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pressures from 1 to 760 Torr [12], In 1977, lYaF-TPI researchers 
demonstrated experimentally that the radius of the plasma channel 
created by an IREB was greater than the beam radius in the pressure 
range of 0.1 to 1 Torr, and the difference between the channel and 
beam radii increased with a decrease in pressure [13]. The same Tonus 
accelerator was used with a 4 cm diameter beam injected into a 9 cm 
diameter drift tube 60 cm long. Beam transport was found to be most 
effective in the pressure range of 0.09 to 0.06 Torr; however, long wave 
hose instabilities were observed in this pressure range with lateral 
beam displacement of the order of one beam diameter (Fig. 22a). 

With further increase in background pressure, beam refocusing 
increased in these experiments because of increased current neutraliza¬ 
tion. In the pressure range for optimum beam transport, a glowing 
cylindrical layer appeared around the periphery of the beam (Fig. 22b). 
This layer was assumed to be a plasma channel carrying part of the 
return current of the beam. The channel was found to narrow down 
with increased pressure; and when the beam developed hose instability, 
the cylindrical symmetry of the layer was not disturbed but the whole 
beam was simply displaced laterally. Measurements revealed that no 
high energy electrons moving parallel to the beam flow were present in 
this layer region. From the analysis of data obtained in Refs. 3, 13, 
and 14, the lYaF-TPI researchers deduced that an IREB can be trans¬ 
ported in large volumes, in large diameter insulating channels, and 



Fig. 22—Photographs of backgroimd gas glow [13] 
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through free space by forming two channels for the return current 
flow—one in the area of the electron beam itself and the other in the 
thin cylindrical layer around the beam. The latter shorted the beam 
current to the accelerator while maintaining beam self-focusing. In 
addition, the plasma layer had a stabilizing effect upon the lateral 
oscillations of the beam improving IREB transport. 

A quickly developing instability was also observed when IREB prop¬ 
agated in background gas at the pressure of about 100 Torr, leading to 
a collapse of the electron beam (Fig. 22c). 

Experimental measurements of the energy spectrum were made for a 
1 MeV, 37 kA IREB injected into a 2 m drift tube filled with air at 
various pressures (Fig. 23) [3,15,16]. Figure 23 shows that the beam 
electrons with less than the maximum energy are located in two 
^oups, one at the front of the beam pulse and the other at the end of 
the pulse. The amplitude of the electron current pulse increased with 
increasing electron energy and peaked at the maximum electron energy 
[3]. This beam energy structure maintained its form over the 2 m 
transport distance in the pressure range of 0.6 to 1 Torr. At 20 Torr, 
the energy spectrum changed considerably after 1 m distance (Fig. 
23b). In this case, both the energy peak and the peak at the end of the 
pulse are decreased considerably. With further increase in background 
pressure to 50 Torr, both the energy peak and the peak at the end of 
the pulse have disappeared completely (see Fig. 23c). Further measure¬ 
ments demonstrated the shortening of the current pulse length at pres¬ 
sures greater than 20 Torr (Fig. 24) [3]. This pulse shortening was 
attributed to current cutoff allowing only the front of the beam pulse 
to pass through the dense gas. 

It was also determined that at long beam transport distances the 
pulse current cutoff appears at lower pressures in the drift tube. This 
current cutoff is attributed to beam instability in dense plasma, which, 
in turn, is created by the ionization of the neutral background gas. 
The lYaF-TPI researchers mentioned in 1975 [3] that experiments 
were being carried out at that time to clarify the mechanism and the 
specific instability characteristics that are responsible for the above 
current cutoff. 

The 1977 experiments [15] also showed the existence of an optimum 
pressure of about 1 Torr when the IREB energy spectrum was least 
perturbed during transport. The maximum energy Em and its width at 
half maximum bE/Em were found to be Em ^ 1 MeV and 
AE/Em ^ 20 percent. Typical oscillograms of the electron beam 
current pulse as obtained at the end of 0.8 and 2 m drift tubes at vari¬ 
ous pressures are shown in Fig. 25. 

Beam transport was found to be most efficient when the beam self- 
fields were fully neutralized along the whole length of the beam path. 
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(b) 


(c) 



p = 20 Torr 


p = 50 Torr 


Observed energies: 

1— 535 to 600 keV 4—805 to 870 keV 

2— 625 to 690 keV 5—895 to 960 keV 

3— 705 to 780 keV 6—985 to 1045 keV 


Fig. 23—Energy spectrum time scan of a 1 MeV maximum energy 
IREB with a 37 kA current pulse at a 1 m transport distance at 
different background pressures [3] 


This occurred at 1 Torr, with the form of the IREB pulse practically 
unchanged (Fig. 25a and b, curve #3). The behavior of Em as a func¬ 
tion of pressure and transport length is shown in Fig. 26. 

Away from 1 Torr, Em was observed to be displaced toward smaller 
energies while the beam energy spread increased. This behavior of Em 
is not attributed to energy loss from gas ionization. 

Below Pcpt, the leading edge electrons were scattered by space-charge 
forces. The nose erosion placed the Em electrons at the lea^ng edge of 
the pulse, and not at the location of maximum current as in the case of 
Popt ~ 1 Torr. The current peak was displaced toward energies 
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(a) 


p = 0,6 Torr 


(b) 


p = 20 Torr 


(c) 


p = 50 Ton- 



Fig. 24—Beam current pulse oscillograph trace for IREB transport in 
air at a distance of 1 meter [3] (reference oscillation 
frequency = 100 MHz) 


lower than thus forming the trailing edge of the pulse during the 
acceleration process. This was the reason for the absence of the high 
electron energy component in the beam spectra and the Em shift 
toward lower energies. Above Popt> the front part of the beam prop¬ 
agated more efficiently. Above 5 Torr (Fig. 25b), beam transport ter¬ 
minated before the Em electrons could appear. These factors lead to 
the same spectra change as in the case for pressures lower than p^pt. 

The Soviet researchers concluded that the observed variation in 
beam energy spectra in beams transported in neutral gas at various 
pressures is due to the different propagation effectiveness of the vari¬ 
ous parts of the beam. This effect is maintained by the transient 










28 


(a) 


(b) 




2 


3 


4 


5 


8 


1 - p = 0.05 Torr 

2 - p = 0.1 Torr 

3 - p = 1 Torr 

4 - p = 5 Torr 


5 - p = 10 Torr 

6 - p = 15 Torr 

7 - p = 20 Torr 

8 - p = 50 Torr 


Fig. 25—IREB current pulses at the end of the drift tube (a) 0.8 m and 
(b) 2 m at various background gas pressures [15]. The vertical line 
denotes the location of electrons with energy 


IREB energy structure determined by the acceleration process. To 
mdntain an unchanged energy character of the beam under transport 
it is necessary to insure conditions that allow full neutralization of the 
beam’s self-fields along the whole path of transport [15]. 

Control of the length of the IREB current pulse was investigated in 
1977 using the principle of beam current cutoff during its transport in 
jmutral gas at pressures higher than Popt [17]. These tests were also 
carried out with the Tonus accelerator at 1 MeV, 40 kA, and a drift 
tube 2 meters long. The current pulse length was controlled by varying 
the pressure within the drift tube. The dependence of the pulse length 
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E^.keV 



1 - L = 0.3 m 

2 - L » 0.8 m 

3- L=1.3m 

4- L*2m 

Fig. 26—Value of the energy at the maximum current peak as a 
function of pressure p and transport length L [15] 


on background pressure and transport length is shown in Fig. 27. 
Pulses shorter than 10 ns can easily be obtained by this method. 

Other experiments at lYaF-TPI made in 1977 on plasma channels 
induced by an IREB propagating through gas and air indicated that the 
beam tracked the plasma channel configuration. These experiments 
were made with the Tonus accelerator with a 1 MeV, 35 kA, 60 ns 
beam injected into a drift space filled with air [14]. The most effective 
IREB transport was achieved at a background air pressure of 0.09 to 
0.6 Torr. At the pressure range of 0.14 to 0.17 Torr long wavelength 
“hose” instabilities were observed with a large lateral displacement of 
the beam. At pressures higher than 0.6 Torr, the beam defocused 
because of an increase in beam current neutralization [14]. 

The above research included Usov’s well-known experiment with a 1 
MeV beam in a large gas-filled drift chamber. In a homogeneous 
transverse magnetic field, the radius of curvature of the beam 
trajectory was determined by the 0.5 MeV nose of the pulse, and not by 
its 1 MeV peak energy, thus indicating that the beam could track the 
plasma channel. 

The researchers at lYaF-TPI stress the importance of the effect of 
collective fields generated by the beam itself and the beam’s interaction 
with the plasma channel. These researchers also found that an IREB 
can be effectively transported across a magnetic field along localized 
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Fig. 27—Pulse length as a function of background pressure and 
transport length, L [17] 


plasma channels having diameters smaller than that of the electron 
beam as well as along thin conducting wires. The effectiveness of 
transport was found to be determined by the degree of magnetic and 
electrostatic neutralization of the beam [14]. 



III. lAE RESEARCH ON HOSE INSTABILITY 


Researchers at the Kurchatov Atomic Energy Institute (lAE), under 
the leadership of L. I. Rudakov, carried out experiments from 1974 to 
1978 using the MS accelerator at 0.4 MeV, 40 kA, 40 ns to investigate 
beam transport and the formation of plasma channels with the help of 
external magnetic fields for beam focusing. The experiments demon¬ 
strated poor beam transport at background air pressures higher than 
0.5 Torr [18,19]. 

Since 1982, lAE researchers have studied the resistive hose instabil¬ 
ity of an IREB beam injected into air and gases [20,21]. Resistive hose 
instability is the result of resonant interaction of magnetic field oscilla¬ 
tions with betatron oscillations of beam particles and is observed as 
macroscopic transverse beam oscillations, which eventually deflect the 
beam completely off axis and onto the wall of the drift tube. 

The formation time of hose instability was found to be considerably 
shorter than the characteristic time of plasma magnetic field decay. 
The distance at which the instability begins to oscillate is of the order 
of several betatron oscillation wave lengths. The time and length of 
the instability formation decreased with increasing pressure from 1 to 
100 Torr. 

Experiments with resistive hose instability were carried out using 
the “Neptun” accelerator at 450 keV, 20 kA, and pulse length of 60 ns. 
The beam electrons passed through a 20 micrometer thick titanium foil 
and were injected into a drift chamber filled with N 2 , H 2 , or air at vari¬ 
ous pressures. The beam transport distance was varied from 3 cm to 3 
m and the pressure within the chamber was varied from 1 to 100 Torr 
[20,21]. 

The measurements involved radial distribution of total beam current 
density at different distances from the point of beam injection at dif¬ 
ferent pressures. Typical scans are shown in Fig. 28 [20]. At 1 Torr, 
the current density distribution maintained its axial symmetry to a dis¬ 
tance of 0.8 m (Fig. 28a). As beam transport distance increased, this 
symmetry was perturbed (Fig. 28b) first at the trailing edge of the 
pulse. In some cases, beam deflection was less than the beam radius at 
a distance of 2 meters, but at 3 meters at the pulse trailing edge, the 
beam was deflected by more than the beam radius. As gas pressure 
increased, the asymmetry of the beam current density profile and beam 
deflection also increased (Fig. 28c). 
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a - p = 1 Torr, L = 26 cm 
b “ p - 1 Torr, L = 85 cm 
c - p = 5 Torr, L = 85 cm 

Fig. 28—Beam current density radial distribution at different times 

after injection [20] 


The measurements [17,18] showed that the distance at which the 
beeun is thrown from the axis into the chamber wall depends upon 
pressure (Fig. 29) and decreases with increasing pressures. 

The data obtained from Faraday cups were in good agreement with 
optical observations of plasma scintillation in the drift chamber 
obtained from electron-optical converter photographs (Fig. 30). 

Figure 30a shows three consecutive photographs of the beam thrown 
off axis onto the chamber wall at 30 Torr. In 10 to 15 ns, the beam 
remained close to the axis of the chamber; in 15 to 20 ns, the beam 
moved radially through about one radius; and in 20 to 25 ns, the beam 
was thrown completely onto the chamber wall. At 60 Torr (Fig. 30b) 
the beam is thrown to the chamber wall in only 11 ns. The charac¬ 
teristic excitation length was 45 cm at 30 Torr and 30 cm at 60 Torr. 
The average transverse displacement velocity was found to increase 
with pressure, changing from 5 x 10® to 2 x 10® cm/s, for a pressure 
change 5 to 60 Torr, respectively [20]. 

The experiments demonstrated that the beam trajectory instability 
develops in a spiral that unwinds toward the end of the drift tube and 
that beam as 5 mimetry is accompanied by beam rotation around its axis. 
The angular velocity of this rotation was 3 x 10® sec”^ and was 
observed in both left and right directions. The instability was found to 
possess a weakly convective character with an oscillation drift velocity 
less than 10^ cm/s. At pressures greater than 10 Torr, the transverse 
oscillation amplitudes were found to be so high that a large part of the 
beam was thrown to the chamber wall. 
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Fig. 29—Distance to off-axis displacement as a function 
of pressure [20] 



a - p = 30 Torr, L = 1.3 m The arrow represents the direction of the 

b - p * 60 Torr, L = 1.3 m beam Injection. The horizontal dotted lines 

represent the position of the chamber wall. 
The vertical dark lines are axial distance 
markers located 20 cm apart 

Fig. 30—Electron-optical converter graphs of axial displacement [20] 
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lAE researchers state that the accuracy of their time measurement 
in the observations of plasma scintillation due to the propagation of 
the beam was within 1 ns. The plasma channel formed by the beam 
was observed to remain for more than 100 ns after the passage of the 
electron beam pulse [20]. 

The current density distribution in the cross sectional plane of the 
transport chamber was measured as a function of time (Fig, 31) [20], 
The figure shows the asymmetry, displacement, and rotation of the 
beam in the cross sectional view. Oscillograms of signals from a beam 
collector placed at the location where the beam was thrown from the 
axis (at L = 85 cm) at p = 15 Torr are shown in Fig. 32 [20]. 


t^ZOns jOns 



Fig. 31—^Beam current density distribution (in A/cm^) as a function of 

time [20] 


The beam propagated symmetrically during the first 15 ns; for the 
next 10 ns the beam disappeared from the collector’s field of view and 
reappeared on one side of the collector’s face. The beam was displaced 
radially with a velocity of 10® cm/s. 

The time of plasma magnetic field decay in these pressure ranges 
was longer than 100 ns for a 3 cm beam radius and plasma conduc- 
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p = 15 Torr, L = 85 cm 

Fig. 32—Oscillograms of the voltage V, diode current I, and beam 
current density as measured with sectioned Faraday cups 
located at different distances from the chamber axis [20] 

tivity a > The beam instability developed in considerably less 

time (t = 1 to 20 ns) than the magnetic field decay time. 

Recent theoretical work at lAE on the resistive hose instability led 
to the conclusion that under conditions of above experiments [20,21], 

^In the investigation of the plasma magnetic field decay time, Soviet researchers do 
not differentiate between the monopole time constant and the dipole time constant 
Tj. The time constant to which the Soviets refer with a value of 100 ns is probably 
which is observed for the plasma current at the tail of a beam current pulse. The value 
of the dipole decay time constant would be 5 or 10 times shorter. The hose instability 
builds up on the time scale of 
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an external axial magnetic field of 6 kG applied to the IREB can 
totally dampen short-wave oscillations causing hose instability. The 
long-wave oscillations do not have time to develop for the meter propa¬ 
gation length conditions in these experiments [22]. 

Theoretical analysis was also applied to beam self-focusing produced 
by two effects in plasma; (1) increase of beam space charge neutraliza¬ 
tion with increasing plasma density and (2) generation of magnetic 
field by thermo-convective currents in inhomogeneous plasma. lAE 
researchers made a calculation to demonstrate the degree of compres¬ 
sion of a beam with parameters of 5 MeV, 60 kA, r ^ 3 cm, in a plasma 
with a density of 10^® cm“^ that extends along the beam direction for a 
distance of L 50 cm. The generation of the magnetic field during 
plasma heating in the beam relaxation process was found to lead to 
beam compression with a compression factor of 10 [23]. 

Other recent theoretical work at lAE included the research on the 
resonant helical instability of an IREB [24] and the study of the gen¬ 
eration of a magnetic field in combination with the filamentation of an 
IREB propagating through a plasma with strong current neutralization 

[25] . These studies showed that for typical experimental parameters 
for IREB passing through gases and plasma (for beam currents from 30 
to LOO kA, = 1 to 2 cm, and the plasma density of about 10^^) helical 
instability is present and fast oscillations are excited causing local 
compression of the beam [25]. 

The expulsion of magnetic field from the electron beam by hot elec¬ 
tron gas whose pressure is greater than the pressure of the magnetic 
field has also been recently investigated [26]. In this theoretical 
analysis, the expulsion process in a high-conductivity medium was 
shown to be accompanied by a rearrangement of the magnetic field in a 
time considerably shorter than skin effect time Tgk {rsk = 4^(rrVc^, 
where c = conductivity, r -= characteristic size of the current flow area) 

[26] . Nonlinear wave solutions have been obtained for a constant 
dispersion velocity, which can be considerably lower than the velocity 
of fast electrons. 



IV. FIAN-IOF RESEARCH—EXPERIMENTAL 
BEAM PROPAGATION STUDIES 


During the early 1970s, experimental work at the Lebedev Physics 
Institute (FIAN) involved general properties during beam transport 
through gases. Since 1981 FIAN, together with the newly formed 
Institute of General Physics (lOF) under the leadership of A. A. 
Rukhadze, has been publishing experimental studies of beam current 
enhancement, thermal effects within the IREB, and their effect on the 
background gas density. 

The first experimental results concerning the perturbation of air 
density by an IREB pulse were published by FIAN-IOF in 1984 [27]. 
The gas dispersion after the initiation of the beam pulse was analyzed 
numerically [28]. In the latest experiment beam parameters were 
measured and gas density change was determined along the beam axis 
[29]. A simple relationship was found between the measured gas den¬ 
sity distribution after the beam pulse and the input energy profile [29]. 

The behavior of gas density after the injection of an IREB pulse was 
observed [27] experimentally using the Terek-IR accelerator [30,31] 
with beam currents of 10 to 15 kA, beam energies of 1.5 MeV, and 
pulse lengths of about 60 ns. The current densities on axis were 
0.3 kA/cm^ and 1 kA/cm^ with the beam radius of 1.5 to 3 cm. The 
beam was injected into atmospheric air through a 50 thick titanium 
foil. 

The gas density radial distribution for different times after beam 
injection is represented in Fig. 33 .[27]. There were considerable differ¬ 
ences in the behavior of background gas with different beam current 
density values. This type of gas behavior could be explained simply by 
the difference in energy input; FIAN-IOF researchers attribute it to the 
energy accumulated in the internal degrees of freedom of the gas and 
to the strong dependence of relaxation time on gas temperature. The 
relaxation time was said to be determined by the input energy of the 
beam. This gas behavior can also be seen in the dependence on time of 
the rarified gas density along the beam axis after beam injection as 
shown in Fig. 34 [27]. 

The dependence of the maximum amplitude of the gas compression 
coordinate upon the time after the beam injection with a current den¬ 
sity along the beam axis of j =1 kA/cm^ is represented in Fig. 35 [27]. 

These experiments determined that the velocity of the compression 
wave outside the volume occupied by the beam itself (r > 3 cm) is close 
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1 - j = 0.3 kA/cm^ 

2 - j = 1 kA/cm^ 

Fig. 34—Gas density along the axis of the beam as a function of time 
after the IREB injection [27] 



t,ns 

Fig. 35—Maximum amplitude of the gas compression coordinate as a 
function of time after the beam injection for j = 1 kA/cm^ [27] 
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to the velocity of sound. The investigation of the density behavior 
after IREB injection can be a valuable source of information about the 
relaxation kinetics of the oscillation-induced conditions. 

The decrease in gas density along the beam trajectory was further 
experimentally investigated [29] using the Terek-IR accelerator with 
beam currents from 4 to 15 kA at beam energies from 1 to 1.2 MeV 
and pulse lengths of about 60 ns. A relationship was determined 
between gas density distribution after the IREB pulse and the energy 
deposition profile. A typical gas density distribution 5p(r) for 
^observed = 34 /iS is shown in Fig. 36a [29]. The radial distribution of 
the beam energy release into the gas is shown in Fig. 36b [29]. 



Solid line - calculation 

Data Points • + - experimental observation 


Fig. 36—a - Gas density radial distribution observed 34 /is after the 
IREB pulse [29] 

b - Radial distribution of the beam energy release into the gas 


The gas density change along the beam axis as a function of time is 
shown in Fig. 37 [29]. The minimum value of gas density near the 
beam axis can be seen from Fig. 38 where its dependence on p/po is in 
good agreement with theoretical calculations. Another issue was the 
total beam current enhancement when the plasma electron current 
flows in the s ame d irection the ffiEB an^ its t^al amplitude c an 
exceed the amplitude of the injected beam. Researchers at Yefremov 
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Solid line - theory 
o - experimental observations 

Fig. 37—Dependence of the gas density along the beam axis 
upon time [29] 



Solid line represents calculated values 

Data points represent experimental observations made in Ref. 27 

Fig. 38—Dependence of the released energy density along the beam 
axis upon the change in gas density at t^bserved >> [^9] 

here: tg ^ where a = effective beam radius 
Q - velocity of sound 
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(NIIEFA) observed the amplification of the beam current in 1978-1979 
[32,33] and they were able to correlate this experimental finding with 
the movement of plasma electrons in the turbulent fields of a plasma- 
beam instability. 

The current enhancement was experimentally determined 
[32,33,34,35] to be due to two-stream instability. However, in 1984, 
FIAN-IOF researchers observed beam current enhancement at pres¬ 
sures close to atmospheric, contrary to what these experiments led 
them to believe [36]. 

Again using the Terek-lR, with 7 kA, 1 MeV, and 60 ns pulse 
length, the beam was injected into air. He, SFg, and a mixture of 
SFe + He -h H 2 . The experiments were carried out both with an axial 
magnetic field {B ^ 2,4 kG) and without a magnetic field. 

The experiments have shown nanosecond spikes with amplitudes 
several times higher than that of the injected beam. This effect was 
observed for all gases at all distances from the point of injection. 

The characteristic length of the spike in the total beam current 
pulse was found to correspond to the time of beam neutralization [30]. 
This was assumed to be the effect of current enhancement at the beam 
front. The spike structure of the pulse was visible with pressure start¬ 
ing from several tens of Torr at distances greater than 1 meter with 
the spike amplitude increasing toward the end of the pulse [36]. The 
application of the axial magnetic field completely obliterated this effect 
even at pressures close to atmospheric. The effect was most prominent 
in helium, but it was also observed in the other gases. 

The effect of total beam current enhancement was attributed to the 
appearance of the spikes in the current pulse. At fairly high pressures 
and at transport distances of greater than 1 meter, at the end of the 
total current pulse a nanosecond peak was observed that was a factor 
of 1.5 greater than the amplitude of the injected beam current. This 
effect of beam current enhancement increased with transport distance 
and with increase in background gas pressure (Fig. 39) [36], 

In He, even at p = 750 Torr, the peaks could not be dampened with 
the axial magnetic field, although other gases did not show this effect. 
There was thus a difference in the conductance of the beam plasma in 
He than in the other gases. FIAN-IOF researchers are at present 
investigating this beam current enhancement theoretically [36]. 

Other experiments performed with Terek-IR included plasma return 
current induced by IREB and concentration and temperature of beam- 
produced plasma, using laser dispersion [37]. At pressures greater than 
3 Torr, plasma concentration increased to 5 x 10^® cm“^ at a tempera¬ 
ture of T -- 1 eV, which demonstrates the presence of collective 
processes that take place in the plasma. For these experiments the 
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Solid line « 0 
Dotted line = 2.4 kG 


Fig. 39—a - Total beam current as a function of distance from 
injection at a pressure of 75 Torr in He [36] 
b - Total beam current as a function of pressure in He at a 
distance of 113 cm from injection 


accelerator was operated with a beam current of 25 kA at 1,5 MeV with 
a 60 ns pulse length. 

The Terek-lR was used in 1984 in a set of experiments on the 
investigation of beam pulse front erosion by neutral gases at pressures 
from 100 to 750 Torr [38], These experiments were performed at 1.5 
MeV, 12 kA, 40 ns using a drift tube 1,2 m long. The gases used were 
nitrogen and a mixture of N 2 and SFq (in the ratio of 6 to 1.) The 
beam transport through the drift tube was carried out with and without 
the application of an axial magnetic field of about 5.8 kG. Erosion of 
the leading front of the beam pulse was observed to increase with beam 
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transport length. At a distance of 60 cm, the pulse width was short¬ 
ened to half its original width. With the application of the axial mag¬ 
netic field this effect of pulse front erosion completely disappeared. 
The erosion of the pulse front was observed over the whole pressure 
range from 100 to 750 Torr, and any quantitative effect of the erosion 
was found not to depend upon gas pressure [38], 

Another set of experiments was performed in 1984 on the injection 
of IREBs into gas at these high pressures but with lower beam current 
densities (75 ^ 50 A/cm^). The Terek-2 accelerator was used for these 
investigations with beam energy of 360 keV and beam currents of 1.4 
kA and 250 A [39]. 

In experiments performed earlier at FI AN using the Terek-lR [30], 
a dependence of the plasma return current upon the beam current as a 
function of background gas pressures was obtained. In subsequent 
experiments [37] the reverse current was found to have a maximum in 
the neighborhood of p - 1 Torr just where a considerable concentra¬ 
tion of plasma electrons was found (Fig. 40), [37]. As pressure 
increased, the role of cumulative ionization (ionization by secondary 
electrons) decreased and direct gas ionization by the IREB became 
responsible for the formation of the plasma. 

The FIAN-IOF researchers stressed the significance of their research 
to optimal IREB transport through gas, as well as to high-power gas 
lasers and plasmo-chemical reactors. The results were also claimed 
applicable to IREB diagnostics in gas, where the measurements of the 
electron beam parameters give no single-valued results or perturb the 
distribution of induced currents or fields. 



+ - He 

• -air 
o-SFg 

Fig. 40—Plasma current to beam current ratio as a function of gas 

pressure [37] 



V, KhFTI EXPERIMENTAL RESEARCH ON 
BEAM PROPAGATION 


KhFTFs experimental research on the propagation of IREBs 
through background gas or plasma can be broken down into two 
separate groups. One group, whose publications extend from 1973 to 
1982 under the leadership of Yu. V, Tkach and Ya. B. Faynberg, has 
investigated IREB propagation through air using a 1 MeV, 50 kA, 30 
ns electron beam with transport distances of up to 3 m. The second 
group, headed by Ya. B. Faynberg with A. K. Berezin as a leading 
author, and publishing from the early 1970s to 1981, used a 2 MeV and 
a 20 MeV linear electron accelerator with low beam currents of 0,5 to 
1.0 A and investigated the interaction of IREB with air and plasma. 
No papers of KhFTI researchers have been found since 1982. 

The early work of the first group, dating back to 1973, dealt with the 
formation of return current, the optimum pressure for transport in air, 
and the measurement of the ionization front velocity in various gases 
in the pressure range from 5 x 10“^ to 10 Torr [40]. In 1980, the 
effect of collisions upon the beam relaxation and of collective and 
elementary processes upon energy transport through gas were studied 
[41]. The later experiments used a 0.7 MeV, 30 to 35 kA electron 
beam, 40 ns pulse length, with transport distances up to 3 m in air, in 
the pressure range from 0.1 to 100 Torr, without the presence of exter¬ 
nal magnetic fields. 

The KhFTI researchers claim to have demonstrated experimentally 
that the relaxation of an IREB in a neutral gas and plasma without 
external magnetic fields occurs mainly through collective processes, 
which, during their formation, are accompanied by the excitation of an 
intensive microwave and X-ray radiation from the plasma, the modula¬ 
tion of the IREB, and the formation of a fast group of particles with 
energies reaching 20 to 30 keV. According to computer simulation, in 
addition to the collective processes, an important role is played by elas¬ 
tic scattering of the beam particles on the plasma atoms and ions, 
which removes particles from the beam channel if the scattering angle 
is greater than a certain critical value. 

The KhFTI researchers devised complex diagnostics for these exper¬ 
iments to obtain accurate measurements. A series of Rogowski loops 
monitored the diode and beam currents at various locations along the 
drift tube, a 25-strip calorimeter measured the total beam energy and 
its radial distribution and allowed control of the beam diameter while 
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making calorimetric measurements, a 12-channel magnetic analyzer 
determined the spectral distribution of the beam electrons, and a 4- 
channel X-ray spectrometer measured bremsstrahlung radiation 
emanating from the plasma heated by the IREB. 

The KhFTI researchers stressed the importance of transporting 
energy to large distances and mentioned that the most important 
characteristic is the beam energy density obtained on target as a func¬ 
tion of background gas. The measured energy carried by the beam over 
different distances as a function of pressure is shown in Fig. 41.. The 
energy carried by the beam as a function of distance for various pres¬ 
sures is shown in Fig, 42.. The condition for maximum energy carried 
by the beam to a certain distance is shown in Fig. 43 as a function of 
pressure [41]. 

From their experimental measurements, the Soviet researchers con¬ 
cluded that in order to obtain the maximum energy carried by the 
beam for short distances ( :< 50 cm) it is necessary to choose the pres¬ 
sure at which the beam is self-focusing at p 1 Torr. To carry the 
beam over long distances, it is necessary to take beam energy dissipa¬ 
tion into account. The optimum gas pressures for long distance beam 
transport were found between 0.2 and 0.6 Torr [41], 

The changes in particle velocity and particle density in beam propa¬ 
gation through gas were measured experimentally. The beam energy 
spread was measured at the output of the accelerator and at various 
locations along the drift tube as a function of gas pressure and was 



= 30 kA, = 1.8 cm, = 5.5 cm; z is path length: 

1- z = 6cm 3-z=1m 

2- z = .5m 4-z = 3m 

Eq = maximum measured beam energy 

Fig. 41—Energy carried by the beam to different distances as a 
function of gas pressure [41] 
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Ijj = 30 kA, Tjj = 1.8 cm, = 5.5 cm 

1 - p = 0.1 Torr 3 - p = 1 Torr 

2 - p = 0.4 Torr 4 - p = 10 Torr 


Fig. 42—Energy carried by the beam as a function of distance for 
various gas pressures [41] 





lb = 30 kA, = 3 cm, rb = 1.8 cm, 

= 5.5 cm 


Fig. 43—Pressure required for maximum energy carried by the beam to 
different distances [41] 
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represented as a time function for two pressures in Fig. 44 [41]. At the 
optimum transport pressure of p = 0.4 Torr, beam energy distribution 
changed only slightly with change in transport distance (see Fig. 44). 
At high pressures, the beam energy peak was shifted to the lower ener¬ 
gies and the low energy particles decreased with increased z, while the 
pulse shortened. 

The beam energy spread was determined as a function of gas pres¬ 
sure for several transport distances (Fig. 45), [41]. The energy spread 
was independent of p and z at short transport distances, but it was 
markedly dependent at the 3 m point. 

Further analysis of their data allowed the researchers to calculate 
the specific beam particle energy as shown in Fig, 46 demonstrating a 
much larger beam energy loss in the pressure region of 0.1 and 1 Torr 
than in the region of 0.3 to 0.6 Torr. 

The researchers noted that at the pressure range of 0.1 to 0.2 Torr, 
the increase in energy loss was due to incomplete space charge neutral¬ 
ization and deceleration of the beam by the space charge field. The 
increased loss in the 1 to 3 Torr range was due to collective processes 
accompanied by microwave radiation. The dependence of microwave 


p = 0.4 Torr 



t, ns 


n, rel. units t, ns 


p-4 Torr 




z = 6 cm 


z = 0.3 m 


z = 3m 


80 ZWJW (rlO 550 E^keV 


Fig. 44—Time plots showing the IREB energy spread at three locations 
along the drift tube for p = 0.4 and 4 Torr [41] 
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n, reL units 



1 - p =s 0.1 Torr 2 - p = 0.4 Torr 

3 - p = 1 Torr 


Fig. 45—Beam energy spectra as a function of background pressure at 
two transport distances f411 


^t/n, reL units 



1-z = 6cm 2-z = 30cm 3--z = 3m 

here is beam energy and n is the number of beam particles 


Fig. 46—Specific energy of the transported electron beam at different 
transport distances as a function of gas pressure [41] 
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radiation intensity upon pressure and transport length is shown in Fig. 
47 for the wavelengths of 3 cm, 8 mm, and 1.8 mm. 

Figure 47 shows that a microwave radiation peak appears at points 
in the pressure range where the beam sustains a considerable particle 
energy loss. The microwave radiation shows a minimum at the point 
of optimum beam transport. As the beam path lengthens, the range of 
pressures for which microwave radiation is observed increases. 



microwave radiation wavelength: 

- • - 4 - X « 3 cm 

——————X = 3 mm 

— — — —-X * 1.8 mm 

Fig. 47—Dependence of the microwave radiation intensity upon gas 
pressure and transport length [41] 


X-ray radiation emanating from the plasma was found to be fairly 
weak in the pressure range corresponding to optimum beam transport 
and to increase at pressures from 0,6 to 3 Torr in the region where the 
maximum modulation of the beam current was observed. The depen¬ 
dence of the X-ray radiation intensity as a function of gas pressure is 
shown in Fig. 48 [41]. 

The angular and radial beam current distribution as a function of 
gas pressure and time were investigated [42] using a beam current of 
30 kA, a beam energy of 0.5 MeV, pulse length of 45 ns, beam radius of 
1.8 cm, and background gas pressure from 5 x 10“^ to 50 Torr. The 
findings are shown in Fig. 49. The initially cylindrical radial electron 
distribution assumed a hollow-beam character at a distance of 40 cm 
from the point of beam injection. At low pressures (less than 2 Torr) 
this happens in a time equivalent to the pulse length; at higher pres- 
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Ip, rel. units 



Fig. 48—Dependence of the soft X-ray radiation (15 to 35 keV) 
emanating from the plasma upon background gas [41] 



a, degrees 


1 - 15 ns, 2 - 30 ns, 3 - 50 ns 
I - 7 X 10*2 Torr II - 2 x lO''* Torr III - 2 Torr 


Fig. 49—Radial distribution of beam current density (a) and IREB 
particle angle distribution (b) in the drift tube for various 
times and pressures [42] 
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sures, the beam possesses a well-defined hollow beam structure during 
the entire pulse. 

Using the 2 MeV and the 20 MeV linear accelerators with a current 
range of 0.5 to 1.0 A, the second group of KhFTI researchers studied 
the interaction of an IREB with plasma. The early experimental 
research of this group (1971-1972) concentrated on IREB interaction 
with low density hydrogen plasma of 10^^ to 10^^ cm~^ using the 2 MeV 
electron linear accelerator at low beam currents of about 1 A 
[43,44,45]- Later experiments (1974 to 1976) involved high density 
plasma of 10^^ to 10^'^ cm"^ [46,47,48]. In experiments carried out in 
1979 with the 2 MeV IREB in air at atmospheric pressure, KhFTI 
researchers demonstrated the formation of a intermediately dense 
plasma (up to 10^^ cm'’^) and also observed an effective collective 
interaction of the electron beam with plasma [49]. The latter was 
apparent from the considerable loss of beam energy as well as from 
microwave and soft X-ray radiation localized in the area of beam 
interaction. 

In 1978-1979 the group published two papers on the LU-40 linear 
accelerator with 20 MeV, 0.5 A, and pulse length of 10 ns [50,51]. 
Only about 1 percent of the beam passed through plasma with a den¬ 
sity of 6 to 8 X 10^® cm“^. At 2 to 3 x 10^® cm“®, about 20 percent of 
the beam could be propagated. However, at a density of 5 x 10^^ cm"® 
only about 15 to 20 percent of the beam would be lost. 

The energy of the X-ray radiation emanating from plasma varied 
from 40 to 90 keV, depending on plasma density, with an intensity of 
lO^jR/cm^-s at 10 cm from the system axis. The Soviet researchers 
claim that the X-ray radiation measurements correlated well with the 
observed beam current losses. 

The interaction zone between the plasma and beam was less than 10 
cm; the beam electrons accelerated by an additional 4 MeV (for an initial 
beam of 20 MeV). The researchers therefore proposed that high- 
frequency fields of at least 400 kV/cm have been set up in the interaction 
region. They have also deduced that the effectiveness of the collective 
interactions of a monoenergetic, low angular divergence IREB increases 
with the plasma density [50,51]. 



VI. NIIEFA EXPERIMENTAL RESEARCH— 
TWO-STREAM INSTABILITY 


NIIEFA’s experiments, under the leadership of L, V. Dubovoy, pro¬ 
vided important information on "reverse current enhancement” (two- 
stream instability) in the propagation of an IREB through air. The 
research was published in 1977 [32,33] with references to earlier 
research dating back to 1973 [52,63] with no new publications since 
1977. 

A 1.5 MeV, 50 kA, 16 ns electron beam from the REP-5 accelerator 
was injected into air and nitrogen at pressures from 10“^ to 760 Torr, 
Most of the investigations were carried out well below the pressure of 1 
Torr where the two-stream instability comes into play. But measure¬ 
ments were also performed at atmospheric pressure in air. The charac¬ 
teristics of the propagating beam were measured at 2 : = 60 cm and 
1.2 m. 

The investigation focused on the case in which the amplitude of the 
current induced in plasma by the IREB pulse was higher than the 
beam current pulse amplitude. This condition was called reverse 
current enhancement (or RCE effect). The RCE effect was present 
only in the case where the electron beam has a large energy spread. 

A typical electron beam current pulse scan is shown in Fig. 50 
demonstrating the beam current total current and reverse 


l,ita 



Fig. 50—^Beam current pulse scan [33] 


53 



54 


current 4(t), where Ir(t) = I^it) - I^it) for a pressure of 0.05 Torr 
[33]. Figure 51 demonstrates the dependence of the basic parameters 
observed in this reverse current process upon the background gas pres¬ 
sure [33]. 

From the experimental results, NIIEFA researchers noted that the 
reverse current can be larger than the beam current 1^ with the 
maximum (as shown for maximum y in Fig. 51) found at a pressure of 
0.03 Torr marked by a maximum in the beam energy spread (AE) and 
a maximum loss in the energy of the beam. 



(a) 


(b) 


(c) 


• -z = 60cm V-z = 1.2m 
a - beam energy E{f )p, here o - beam radius r^{f)p 
b-Y(p) = (l/l,)(p) 

c - AE(p) change in energy of the beam spectra 

Fig. 51—Dependence of system parameters upon background 
gas pressure [33] 







VII. NIITP EXPERIMENTAL RESEARCH—IREB 
INJECTION INTO ATMOSPHERIC AIR 


From 1977 to 1981 the NIITP carried on a seemingly minor effort 
(represented by two published papers and four authors) under the 
leadership of V. M. lyevlev. A 500 keV CW electron beam was injected 
into air at atmospheric pressure through a differential pumping system 
instead of a titanium foil window [54,55]. Special attention was given 
to the development of the differential pumping system utilizing a series 
of apertures with various pumping schemes incorporated between each 
aperture in order to reduce the gas load and provide a smooth transi¬ 
tion between the atmospheric air and the high vacuum required in the 
accelerator section. 

The maximum propagation distance of 200 keV electrons through 
atmospheric air was only about 30 cm. However, as the air was heated 
up by the beam, a lower density channel formed, extending the propa¬ 
gation distance. No further information was published on the details 
of beam propagation through atmospheric air. 

Stated practical applications of the research were beam welding, cut¬ 
ting, and hole boring of solid materials. 
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VIII. VEI RESEARCH—THEORY AND 
EXPERIMENT 


VEI researchers have recently performed extensive research on 
IREB transport through low-pressure background gases (10”^ to 10”® 
Torr), as described in a previous report [1], In high gas pressure, 
IREB propagation studies that began in 1979 are still in progress. VEI 
researchers, under the leadership of A. V. Zharinov and M. A. Vlasov, 
concentrated on the study of the propagation windows for IREB trans¬ 
port, the determination of the upper limit of gas pressure to maintain a 
plasma-beam discharge (plasma formed by IREB injected into air), 
injection methods of an IREB into the atmosphere using differential 
pumping, the formation of large microwave signals inherent in the 
beam, and the formation of positively charged electron beams. The 
latter of these seems to be of interest because the researchers claim 
that the beam does not become subjected to any irregular instabilities 
and is axially self-centering. 

The upper limit of the gas pressure for a plasma-beam discharge was 
experimentally determined in [56] at VEI using a 300 keV, 20 A, 400 
fiSy 0.6 cm diameter electron beam and injecting the beam into air 
through a differential pumping system. The pressure range used was 
0.1 to 300 Torr. It was determined that in the pressure range 50 to 300 
Torr the electron beam propagation through the dense air was stable 
and that the dynamics of the beam radius change along the length and 
time of propagation can be well approximated by the gas collision 
theory. 

Researchers at VEI have investigated the decrease in the number of 
neutral gas atoms in the path of the plasma channel due to the effect 
of the hydrodynamic expulsion of the channel hot gas [57,56]. Here 
the interaction of the IREB with the background gas has been investi¬ 
gated showing the basic mechanism to be the dissipation of energy due 
to electron collisions with the background gas particles whereby plasma 
is formed and the gas is heated up. 

The gas temperature was estimated from the following expression 

T ~ (p„/6)®-2 

where Pa = the initial gas density 
Ih = beam current. 


56 



57 


For pressures of 10 to 40 Torr, the value of T was found to be 0.4 to 
0.5 eV for the observed beam parameters; the concentration of neutrals 
inside the channel was calculated to be 2.5 to 8 x 10^® cm"^ [56]. 

VEI researchers carried out experiments in 1984 [58] on the forma¬ 
tion of a beam-plasma discharge in the area where the electron beam 
was injected into the atmosphere through a differential pumping sys¬ 
tem. The energy of the beam was 300 keV, beam current 27 A, pulse 
length 500 and diameter 0.6 cm. The beam transport was stable 
through the differentially pumped beam injection area where the 
plasma density had an extended gradient and the gas pressure changed 
from 20 to 400 Torr. 

The interaction of the electron beam with the background gas in 
this area led to the heating of the gas and the gas radial expansion, 
which generated a cylindrical shock wave. The latter was reflected off 
the walls, collapsed onto the axis, and radially expanded. Because of 
the gas heating, a gas plug was formed in the channel and obstructed 
the gas flow. It was demonstrated that because of the conditions set 
up by the gas density gradient along the beam axis, a plasma-beam 
discharge was created with microwave radiation in the 3 cm wavelength 
region. 

VEI researchers in 1984 also described a differential pumping system 
they have designed for the injection of an electron beam from a pres¬ 
sure of 10“® Torr to a region under pressure of 10 to 100 Torr [59]. 

An interesting finding of the recent VEI research was the formation 
of overall positively charged beams of electrons [60,61]. These beams 
had an overabundance of positive charge in the neighborhood of 10® to 
10^ V and could be varied as required by changing either the pressure 
or the type of background gas or metal vapor. The critical current was 
shown to increase several hundred fold (over 10 A) using gridded con¬ 
ducting tubes for propagation. The electric fields in this case coidd be 
higher than in the quasi-neutral case (when Ui = ne, r£ = Xg) by factors 
of 100 to 1000. The main advantage of positively charged electron 
beams is that they are not subject to irregular instabilities and have no 
beam current limits because of these instabilities and are also axially 
self-centering, which may allow them to be easily transported along 
various configurations of metal beam pipes [60]. The formation of this 
type of electron beam is explained by VEI researchers using the 
apparatus shown in Fig. 52 [60]. The electron beam flows with a beam 
ctirrent /t, velocity u = /3c along the axis of a tube of length L and 
radius R as shown in the figure. A potential Vt is applied to the tube, 
which is negative with respect to other electrodes that contact the 
plasma created by the beam. All the secondary electrons created by 
the beam inside the tube have to leave by the ends, and all the ions 
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Fig. 52—Apparatus for the formation of positively charged electron 
beams showing the voltage radial distribution [60] 

c 


must flow to the walls imder the action of the field where they 
recombine at the wall surface. For the secondary electrons to leave the 
tube, the voltage on the tube must be ^ greater than the voltage 
between the beam boundary and the inside of the tube <pR (as shown in 
the voltage distribution plot in Fig. 52). The condition for the electron 
beam to pass through the tube must be met: Here the kinetic energy 
of the beam ipQ must be larger than The calculations in Ref. 

60 demonstrate that this sort of overcompensated electron beam can 
indeed be formed and could be of future interest because of its 
aforementioned advantages. 

VEI also considered the practical aspects of the electron beam for 
processing of materials [57,62], involving the injection of 0.5 MeV 
beam into air at atmospheric pressure. It was found that electron 
beam propagation is very weakly dependent upon beam current 
(L - but strongly dependent upon beam energy (L - Es^). Beam 
energies of the order of 0.5 MeV were recommended for processing of 
materials. Even at this energy, the propagation distance was about 0.6 
cm. The researchers also mention that energy transport by means of 
an electron beam, while taking into account scattering and gas heating, 
still has not been investigated sufficiently [57]. 





IX. MFTI RESEARCH ON BEAM INSTABILITY 


Researchers of the Moscow Physico-Technical Institute (MFTI), 
under the leadership of K. V, Khodatayev, carried out investigations 
from 1973 to the present on the theory of the resistive hose instability 
and the development of beam oscillations and noise. 

A nonlinear solution to the sausage mode (azimuthal quantum 
number m = 0) of the resistive instability was determined in 1984 [63] 
for an IREB passing through a conducting plasma. It was determined 
that the sausage mode of the instability actually does develop. In the 
case of a Bennett distribution function, the sausage mode of the resis¬ 
tive instability has been observed to be dampened. 

Using different values of exitation frequencies, these calculations 
showed that phase mixing leads to complete damping of the sausage 
mode of the resistive instability in the case of a hot beam with a Ben¬ 
nett density distribution, According to MFTI researchers, this conclu¬ 
sion holds true in conjunction with the experimental results of Refs. 64 
and 65 in which the hose instability with mode m = 1 was observed, 
but where there was no observation of the sausage mode of the resis¬ 
tive instability. 

MFTI researchers have looked at the energy of beam noise as a 
function of time. They have seen this beam noise mount rapidly, reach 
a saturation, and then slowly decrease with time [66]. They have 
determined that the beam distribution gradually stops varying, the high 
frequency oscillations become dampened, and these oscillations transfer 
their energy to the thermal movement of the beam particles. Here the 
beam emittance, which is the characteristic of the thermal energy of 
the beam, begins to increase. 

Equations of nonequilibrium states of a hot ribbon electron beam 
moving through a stationary ion background indicate that the beam 
relaxes to an equilibrium condition. That process was broken down 
into two cases. In the first case, the transformation took place in con¬ 
junction with the collisionless theory in which the phase mixing of 
electron trajectories led to a large modulation of the beam distribution 
function. The calculations showed that this quenches the oscillations 
of the macroscopic beam characteristics. In the second case, reversible 
processes smooth out the inhomogeneities of the distribution function 
in phase space, and thus produce beam relaxation. 

The characteristic time for beam relaxation was determined to be 
the electron oscillation period. The rate of oscillation damping of the 
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macroscopic beam characteristics depends upon the type of initial exci¬ 
tation of the beam distribution. Small amplitude oscillations, which 
represent a fairly small area of the phase space, are damped slower 
than the large excitations, which possess a larger phase space [66]. 



X. lEM RESEARCH 


Researchers at the Institute of Experimental Meteorology (lEM) in 
Obninsk were working on the transport of IREBs through neutral 
gases in 1980-1981 [67,68]. In their theoretical studies, they have 
analyzed a condition of self-channeling of the electron beam where the 
action of the beam’s self-magnetic field is capable of considerably 
decreasing the beam’s spreading due to multiple scattering. In these 
investigations, an example was chosen with a 20 MeV, 34 kA electron 
beam transported in air with a density of 1.3 x 10”"^ g/cm^ with an ini¬ 
tial beam cross section of 100 cm^. It was determined that the cross 
sectional area at 100 meters increased a factor of 4.5. When the beam 
current was taken at low values (/^ = 100 A), however, this value 
increased a factor of 3 for all other conditions unchanged [67]. 

These researchers reached the following conclusions [67]: 

1. The condition of self-channeling is analogous to the equilib¬ 
rium case for the pinch effect where the transverse energy of 
the beam particles acts as the temperature. 

2. The effect of the beam’s self-magnetic field can considerably 
decrease beam spreading caused by multiple scattering. 

3. Even with self-channeling, the beam spreads; initially this fol¬ 
lows an exponential law. 

These observations, however, have been well documented by 
Western sources in the early 1960s and 1970s. 
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XI. CONCLUSIONS 


Soviet research on the propagation of intense relativistic electron 
beams appears to be broad and systematic, and it has gathered momen¬ 
tum during the past five years. One obvious application of such 
research is in the area of inertial fusion reactors driven by electron 
beams. Another is to resolve a problem in the quest for the charged- 
particle beam weapon—the feasibility of efficiently transporting such 
beams through atmospheric air. The data presented here and drawn 
from open Soviet publications are pertinent to some basic issues of 
beam weapon research. However, they cannot be directly attributed to 
such research, primarily because of the low beam energies (<1.5 MeV), 
short propagation distances (< 3 m), and low pressures used in the 
published Soviet experiments. 

An indirect hint of a possible military involvement is the disappear¬ 
ance from published journals of five key experimental teams during the 
period 1978 to 1981, including a lYaF-TPI team lead by Yu. P. Usov 
and two teams at KhFTI led by Ya. B. Faynberg, Yu. V. Tkach, and A. 
K. Berezin (see Fig. 1). If the disappearance of these teams was due to 
security classification, then the total Soviet effort in this research area 
may be much larger than appears in the published literature. Even in 
the absence of publications from the five experimental research teams, 
the overall number of Soviet published papers has steeply increased 
since 1983, and the number of authors has shown a similar increase 
since 1980 (see Fig. 2). 

Five new Soviet research efforts, initiated in the late 1970s or early 
1980s, are currently active in IREB propagation studies. Of these, 
FIAN-IOF, led by A. A. Rukhadze, demonstrates the largest experi¬ 
mental effort; the largest theoretical effort is being carried out at 
lYaF-TPI under the leadership of V. P. Grigor’yev. 

Because the research being carried out by the new research teams 
was basically the same as that of the five old teams, two hypotheses 
are offered: 

1. The entire research program has been transferred to the new 
teams, while the old team membership has been disbanded. 

2. The new teams assumed only the part of the program suitable 
for open publication, and the old teams commenced a new, 
advanced phase of research subject to classification. 
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The second hypothesis is more probable, because of the total absence 
of the previous authors in the literature of the 1980s. It is unlikely 
that disbanded team members would altogether cease contributing to 
scientific research. 

The specific research issues pursued in the published literature by 
the former and the new Soviet research teams on IREB propagation 
through air for the past 15 years are as follows: 

• The formation of a plasma channel created by an IREB propa¬ 
gating through background air and the investigation of the effect 
of beam parameters upon the plasma channel parameters (and 
vice versa). 

• Determination of the background air pressure for the optimum 
transport of an IREB, Two optimum pressure ranges were deter¬ 
mined: 

— The ‘Hon focused regime”—0.06 to 0,09 Torr. 

— The “low pressure window”—1 Torr. 

• The current enhancement effect, whereby the current induced in 
the plasma by the IREB is larger than the original beam current. 

• The effect of the resistive hose instability upon IREB propaga¬ 
tion. 

These key issues of Soviet research match the main interests of U.S. 
research on IREB propagation through high pressure air. 

Soviet researchers found that an IREB can be transported in large 
volume, large diameter insulating channels, and through free space by 
forming two channels for the return current flow: one in the beam and 
the other in a thin cylindrical layer around the beam. The latter 
shorts the beam current to the accelerator while maintaining beam 
self-focusing; and the plasma layer has a stabilizing effect upon the 
lateral oscillations of the beam, thus improving IREB transport. This 
layer was assumed to be a plasma channel carrying part of the return 
current of the beam. When the beam developed hose instability, the 
cylindrical symmetry of the layer was not disturbed, but the entire 
beam was observed to be displaced laterally. Measurements revealed 
that no high energy electrons moving parallel to the beam flow were 
present in this layer region. It was determined that the plasma chan¬ 
nel was produced by the ionization of the background gas by electrons 
of the IREB leading front, and the beam trajectory was determined by 
the movement of the low energy electrons and not by the main part of 
the beam. The plasma channel remained intact longer than 100 ns. 
The beam was also found capable of tracking the plasma channel. 

In their investigation of “current enhancement,” Soviet researchers 
discovered that most of the plasma current density Jp flowed within 
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the confines of the electron beam in a time about 100 ns, which is con¬ 
siderably longer than the beam pulse. The return current enhance¬ 
ment (RCE) effect was found to be present only in the case where the 
IREB had a large energy spread. The limitations of the ‘‘two-stream 
instability*’ at the low pressure of IREB propagation has been reported 
only by NIIEFA, and the earlier Soviet reports have mentioned that all 
current enhancement is due to the two-stream instability. However, in 
1984, FIAN-IOF researchers observed a current enhancement condition 
that takes place at pressures close to atmospheric when none of the 
earlier studied causes of current enhancement could be the reason for 
the effect. This anomalous current enhancement was found to increase 
with pressure and propagation length and was in the direction of the 
IREB. The Soviet researchers claim that they are currently making 
further studies of this effect. 

In the Soviet open literature, the only mention of laser guidance 
(the use of a laser to create a low-pressure channel along whose axis an 
IREB can propagate) was made in 1974 at FIAN-IOF by G. A. 
Askar’yan [69]. Here a preliminary experiment demonstrated the for¬ 
mation by a laser of a low-density channel in air. This channel had a 
finite lifetime of several hundreds of microseconds with a density 10 to 
20 times lower inside the channel. A 50 keV electron beam pulse was 
passed, through this channel successfully, thus demonstrating laser 
guidance. Other FIAN-IOF papers published by G. A. Askar’yan and 
his co-workers have made statements on the use of laser guidance for 
the successful transport of charged particle beams. There are no recent 
publications on laser guidance, and it is believed that this subject is 
censored by the Soviet authorities. 

Most of the Soviet reports that have been considered here do not 
specifically mention a single application of their research. However, 
many papers generally indicate the importance of optimum propagation 
conditions for high energy, high current beams, and their transport to a 
maximum allowable distance. In a few cases (VEI, NIITP, MFTI), 
where the injection of the beam into atmospheric pressure air was con¬ 
sidered, specific applications for beam welding, beam cutting, and 
specimen hole boring were mentioned along with the general objective 
of achieving maximum propagation length in air. 



Appendix 


SOVIET RESEARCH TEAMS INVOLVED IN 
IREB PROPAGATION 


Soviet research teams involved in IREB propagation in air at pres¬ 
sures from 0.1 to 760 Torr are listed below by institute affiliation. 
Group leaders are indicated in italics. 

lYaF-TPI: Institute of Nuclear Physics at Tomsk Polytechnic Institute, Tomsk 
Team 1: Theory and Computer Simulation 
V. P. Grigor'yev 

A. G. Potashev N. S. Shulayev A. V. Zakharov 

A. 1. Koryakin A. S. Popov A. N. Didenko 

G. P. Isayev 

Team 2: Experimental 
Yu. P. Usov 

A. N. Didenko V. A. Tuzov A. I. Ryabchikov 

Ya. Ye, Krasik A. V. Petrov G. 1. Kotlyarevskiy 

lAE: Kurchatov Institute of Atomic Energy, Moscow 
L. L Rudakov 

V. P. Smirnov V. V. Gorev S. V. Zakharov 

L. Ye. Aranchuk V. D. Vikharev S. F. Zhandarov 

V. D. Korolev L. I. Urutskoyev S. F. Grigor'yev 

A. V. Gordeyev V. V. Zazhivikhin V. Yu. Shuvayev 

FLAN: Lebedev Physics Institute/IOF: Institute of General Physics, Moscow 
A. A, Rukhadze 

A. A. Savin G. P. Mkheidze A. L. Ipatov 

G. A. Askar’yan B. M. Manzon R. R. Kikvidze 

I. M. Minayev A. G. Shkvarunets V. Yu. Shafer 

S. N. Kabanov A. A. Korolev A. 1. Kuz’min 

V. Ye. Kul’beda V. Yu. Nosach 

KhFTI: Physical-Technical Institute, Khar’kov 
Team 1 

Ya. B. Faynberg 
Yu. V. Tkach 

1.1. Magda G. V. Skachek S. S. Pushkarev 
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V. A. Bondarenko 
S. L Naysteter 
P. M. Lebedev 

Team 2 

Ya. B. Faynberg 

A, K Berezin 
V. A. Kiselev 

B. G. Safronov 

N. M. Zemlyanskiy 

L. I. Bolotin 
V. A. Buts 


L P. Panchenko 
I. N, Onishchenko 


V. P. Zeydlits 
G. L. Fursov 
1. K. NikoFskiy 
S. S. Krivulya 
V. I. Kurilko 


B. D, Panasenko 
V. I. Shevchenko 


1. A. Grishayev 
L N. Onishchenko 
V. L Mirnyy 
A. M, Yegorov 
A. P. Tolstoluzhskiy 


NIIEFA: Yefremov Institute of Electrophysical Equipment, Leningrad 

L. V, Dubovoy 

T. S. Gosteva G. R, Zablotskaya B. A, Ivanov 

V. I. Chernobrovin S. A. Kolyubakin 

NIITP: Scientific Research Center of Thermal Processes, Moscow 
V. M. lyevlev 

A. S. Koroteyev V. V. Koba I. G. Kulakov 

VEI: All-Union Electric Engineering Institute, Moscow 

A. V. Zharinov 

M. A, Vlasov 

S. V. Nikonov S. L Vybornov V, A. Safonov 

A. V. Rykhlov V, P. Popovich Yu. A. Yel’chin 

V. A, Artemov L. G. Dubas I. P. Denisova 

A, V. Smirnov 


MFTI: Moscow Physico-Technical Institute, Moscow 
K. V. Khodatayev 

V. G. Leyman A. A. Ivanov K, G. Gureyev 

V, 0, Zolotarev S. D. Stolbetsov N, Ye, Rozanov 

N. S. Putvinskaya 0, B. Ovsyanikova F. F, Kamenets 

lEM: Institute of Experimental Meteorology, Obninsk 

A. P, Budnik 
P. N. Svirkunov 
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